Entropic competition between knots and slip-links 
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Using canonical Monte Carlo simulations, we introduce a new numerical procedure for comparing 
the entropic exponents of polymers with different constraints and/or topologies. Setting up compe- 
titions between polymer segments which can exchange monomers according to their free energies, 
we obtain the universal exponents of partition functions, independently of any knowledge of the 
non- universal part. The method is successfully tested for closed polymer loops decorated with slid- 
ing rings. We also investigate the entropic exponents of loops with a fixed knot type, in which case 
we are limited by strong finite-size effects. 



I. INTRODUCTION AND SUMMARY 

A long flexible polymer chain can assume an enormous 
number of configurations; thus, the conformations of 
polymers are best described statistically. Despite the mi- 
croscopic distinctions which exist among different poly- 
mers, at large length scales many generic properties are 
independent of the details of specific polymer structure 
and thus are "universal." yj The understanding and char- 
acterization of this universal behavior has led to great 
advances in the statistical treatment of polymers. Poly- 
mers in good solvents have been modeled by self-avoiding 
walks (SAWs), which are in turn mapped to a magnetic 
system at its critical point pj. Renormalization group 
techniques then allow a variety of polymer properties 
to be calculated analytically 0]. On the experimen- 
tal front, advances in manipulation and imaging of single 
molecules have provided new impetus for studies of poly- 
mers With these micromanipulation techniques, it is 
now possible to explore a wide variety of physical proper- 
ties of polymers, and thus to test models and theoretical 
predictions at the level of a single chain. 

In many circumstances, the number of configurations 
of a polymer grows with its length £, as 

w (£)~Aii e r c [i + Be- A + ■■■]. (i) 

In the above equation, the connectivity constant fi, is a 
non-universal quantity which depends on the microscopic 
features of the chain, but (like a free energy density) re- 
mains the same for different boundary and topological 
constraints By contrast, c is a universal exponent, 
which is independent of the microscopic characteristics of 
the polymer, but which does depend on the dimensional- 
ity, and on global boundary and topological constraints. 
For example, in the case of a closed loop, c = dv, where d 
is the dimensionality of space, and v is the exponent re- 
lating the typical spatial extent of a polymer to its length 
by R ~ £ v . (In d = 3, for polymers with self-avoiding 
interactions v w 0.588:1, 4].) In Eq. (JTJ, we have antic- 
ipated subleading corrections to the leading asymptotic 



behavior at large £, indicated by so-called corrections to 
scaling in the square brackets. Renormalization group 
calculations suggest that the exponent A is universal, 
while the amplitude B is case-specific [2 • 

Numerical determination of a power-law correction to 
the leading exponential behavior is rather cumbersome. 
It is usually accomplished by examining a generating 
function f(z) — ^ t £ 9 w(£)z l , with a conveniently chosen 
value of the free parameter g. The function f(z) is usu- 
ally singular at z — l//i [f| , and the details of this singu- 
larity determine the exponent c appearing in Eq. ^ This 
can be done either by examining a finite amount of terms 
in series expansions, using w{£) extracted from exact 
enumeration studies, or by performing grand-canonical 
Monte Carlo(MC) studies of the polymer with weight 
determined by £ 9 z l for several values of z, as was done 
inRefs. H0. 

Here, we propose a method for calculating the uni- 
versal exponent c in the number of configurations, with- 
out the need to account for the leading non-universal 
exponential growth fi . Our procedure is best suited for 
the evaluation of the exponent c when the polymer is 
constrained in various manners, a situation that occurs 
physically when a polymer can move from one regime to 
another. Using canonical Monte Carlo simulations, we 
directly compare the relative number of configurations by 
setting up an ensemble in which two segments of a poly- 
mer can exchange monomers. A pair of monomers be- 
longing to two different segments does not interact, while 
the monomers in the same segment repel each other. 
Each segment is subject to its own set of global con- 
straints. The total number of monomers is fixed to L, 
while in any given configuration the two segments have £ 
and L — I monomers, respectively. According to Eq. IjTjI. 
in the asymptotic limit, the number of configurations of 
such a system is 

w(£,L) = A£- Cl (L-£)- C2 , (2) 

where c\ and C2 are the exponents characterizing each 
segment; the prefactor A is independent of £. Assuming 
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FIG. 1: Examples of possible applications of "entropic com- 
petition:" (The universal power-law exponents depend on the 
boundary condition of the system.) (a) is appropriate for cal- 
culation of the exponent 7 which describes open linear poly- 
mers. Note that the two segments of the chain do not interact 
with each other, (b) is good for calculation of the exponent 
71, (c) for a wedge, (d) for exponent 711, and (e) for a chain 
attached at one end to the interior surface of a sphere. 



that the asymptotic limit has been reached, fitting the 
histogram of t should provide the exponents c\ and c^. 

This method, which we call "entropic competition," 
can in principle be applied to several physical circum- 
stances, some of which are depicted in Fig. ^ For ex- 
ample, one can extract the entropic exponents of linear 
chains by threading them through a hole, and making the 
first and last monomers of the chain bigger than the size 
of the hole such that the chain can diffuse back and forth 
without wandering away (Fig. ^1). A key point in this 
set-up is that the two parts of the chain are not allowed 
to interact with each other. By examining the probabil- 
ity of having a segment of size £, which is proportional 
to £ 7_1 (L — £) 7_1 , one can extract the entropic exponent 
of linear chains, known as 7 ||. If, instead of thread- 
ing a chain through a hole, we divide the chain into two 
parts by a rigid wall and let the two parts of the chain 
exchange monomers at the wall (a process called translo- 
cation), we can extract the entropic exponent which is 
called 71 [j| (Fig. [TJ5). Other potential applications are 
to the calculation of entropic exponents of linear poly- 
mers in the presence of a wedge, or inside a spherical shell 
(Figs.^c, e)). The latter case is particularly relevant to 
the translocation of DNA through spherical capsids and 
has been the subject of intense research^. "Entropic 
competition" is especially appropriate for calculation of 
the entropic exponent 711, related to the polymers with 
two ends anchored to a surface In this case, by fixing 
one monomer of a ring on a wall and letting the two seg- 



ments of the ring on the two sides of the wall exchange 
monomers with each other at any point on the wall, we 
can extract the entropic exponents of chains with two 
ends restricted to move on a surface (Fig. QJ) . The latter 
is of renewed interest due to potential relevance to a new 
class of physical gels obtained from triblock copolymers 
anchored to fluid membranes at each end 

El- 

The calculation of entropic exponents using "entropic 
competition" is not solely restricted to open polymers. 
One can divide a ring into two loops with the use of a 
'hole' connecting two separate spaces, and a monomer 
which is fixed in position, as illustrated in Fig. [21 In 
Sec.|n]we employ this system to calculate the probabil- 
ity distribution of two non-interacting self-avoiding (SA) 
loops of lengths I and L — £ , which can freely exchange 
monomers. If the two polymer segments have the same 
topology and boundary conditions, as in the aforemen- 
tioned case, then the exponents c\ and C2 of Eq. (0) are 
equal and one can directly find the entropic exponents by 
fitting the probability distributions resulting from simu- 
lations to Eq. 0. 

The entropic exponent can in fact be modified in a sys- 
tematic manner by decorating polymer loops with slid- 
ing rings. We investigate the effect of sliding rings on 
the number of configurations of closed chains in Sec. IIHI 
We assume that the size of a ring is the same as that of 
one monomer. In this situation, the presence of a ring 
increases the number of configurations of the loop by a 
factor of I (the number of monomers on the loop where 
the ring can slide). Once again we examine the accu- 
racy of "entropic competition" by comparing the results 
of the simulations with analytical expectations, and find 
an excellent agreement between them. 

Introduction of the "entropic competition" method 
was motivated by its potential application to knots. The 
influence of a knot on the entropy of polymers has been 
the subject of interest for some time 0, 0, El E3- In 
Sec . 11 VI we study knotted polymers. We overview the cur- 
rent hypothesis regarding the number of configurations of 
a knotted polymer. We numerically study several related 
models and observe that the presence of knots changes 
the number of configurations. Furthermore, we numeri- 
cally calculate the probability distribution of two loops 
with one trefoil knot on each side and find that it is com- 
pletely different from that of two loops with one ring on 
each side. This may well be an indication of the impor- 
tance of finite-size effects in knots, as the two cases are 
conjectured to be similar in the limit of very long chains. 

Finally, in Sec. E] we observe the "competition" be- 
tween a knot and rings by placing a knot on one loop 
and rings on the other loop. In our simulations, the knot 
pulls the chain much harder than a single ring, possibly 
due to finite-size effects. For polymers in the range of 
a few hundred monomers, roughly six rings are needed 
to "compete" with a simple trefoil knot. We attempt to 
quantify this finite knot-size effect. 

In the absence of an exact statistical treatment of 
knotted polymers, several researchers have employed 
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FIG. 2: A schematic depiction of two loops in competition. 
Sides (a) and (b) do not interact with each other. The position 
of one bead is always fixed (the one with X inside). The oval 
ring shows the position of the 'hole' separating the spaces in 
which the segments exist. 



'slip-links' for gaining analytical insights on topological 
constraints 1 12l Il3j . Slip- links can be envisioned as belt- 
buckles which force two points of the chain to be close to 
each other. Para-knots are collections of such slip-links, 
and have been proposed as means of designing entropy- 
driven functional molecules which are linked to each other 
mechanically rather than chemically flij. In Sec. we 
also employ a para-knot model to explore the effects of a 
sub-leading scale on entropic competition. 



II. METHOD 

We start by calculating the entropic exponents of a 
closed loop. Our model polymer consists of L hard 
spheres of diameter D connected into a chain by tethers 
which have no additional energy cost 0] , but restrict the 
distance between connected spheres to be smaller than 
1.2D; this prevents the chain from crossing itself. Fig- 
ure |21 provides a schematic depiction of the simulation: 
Solid circles represent monomers of the chain; the po- 
sition of one monomer (labeled with an X) is fixed in 
space, and the chain is forced to pass through a hole 
which is depicted by an oval ring in Fig. The effect of 
the hole combined with the fixed monomer is to create 
two loops of lengths I and L — £. The two loops do not 
interact with each other: as soon as a monomer passes 
through the hole, it ceases to interact with the monomers 
of its previous side and starts to interact with the ones in 
the new side. In order to preserve the topology of each 
loop, we allow the monomers immediately adjacent to the 
hole on either side to interact with each other, and also 
allow the fixed monomer to interact with all the other 
monomers. It is also necessary to keep the hole narrow 
enough not to let knots pass through it (for later appli- 
cations). The hole creates a barrier which slows down 
the process of transferring the chain from one side to the 
other one. In fact, this entropic barrier is the major ob- 
stacle to the computation of the probability distribution 
of long chains (L > 300) within a reasonable amount of 
CPU time. Obviously, this barrier depends on the en- 
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FIG. 3: Probability distribution for two loops exchanging 
monomers. The dots are the results of the simulation for 10 9 
MC steps, with L — 100. The graphs are normalized such 
that the integrated weight is equal to one. 



tropic exponents, and is less significant for a hole on a 
rigid wall. In calculating the exponent 711 we remove 
the entropic barrier (hole) completely; the two segments 
on the two sides of the wall can exchange monomers at 
any point on the wall (Fig. QJ;). In this case, we let the 
monomers pass through the wall only based on their or- 
der along the chain, i.e., monomer n can pass through 
the wall, only if monomer n — 1 has already moved to the 
other side of the wall. 

The number of configurations of the system depicted 
in Fig. |2 is the product of number of configurations of 
each loop, and thus scales as 



w 00 (£, L- t) = A^ L r vd x(L-t) 



(3) 



where I is the number of monomers in one of the loops. 
(The absence of a knot or ring on each loop is denoted 
by the index pair 00.) Note that the n does not affect 
the ^-dependence, and shall in fact be ignored henceforth. 
Although the passage through a worm-hole may appear 
unphysical in this context, our simulations produce the 
correct results for the two non-interacting self-avoiding 
loops which exchange monomers. Figure |3] shows the 
probability distribution for the segment length I. The 
dots are the result of the simulation for 10 9 MC steps fl6) . 
The solid line corresponds to the Eq. with v = 0.588 
at d = 3. The graphs in this figure and all the other 
figures are normalized such that the integrated weight is 
equal to one. The good match of simulation with Eq. © 
confirms that our method for computing the probability 
distribution is accurate. Similar curves were produced 
through simulations for chains with different sizes L = 
50, L = 150, and L = 200; all follow Eq. ©. 



III. SLIDING RINGS 

To further examine the validity of our method, in this 
section we present results of simulations including sliding 
rings. For each configuration of a loop of length £, a 
sliding ring can occupy £ different positions, and thus 
the probability distribution of such a loop scales as, wi ~ 
£wq ~ i/£~ wd+1 . In simulations, the passage of a chain 
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through the hole is stopped if a ring is placed exactly on 
the monomer entering the pore. Instead of simulating 
a "real" sliding ring, we calculate the probability that 
a ring might be sitting adjacent to the hole and then 
prevent the chain from passing through the hole with 
this probability. If we have one ring, the probability of 
the ring to be exactly on the monomer entering the hole 
is l/£. In case of many rings, this probability goes up. 

We start by placing a ring on one of the loops of Fig.|3 
and then we measure the histogram of lengths. Figure 0] 
shows the probability distribution of the loop side with 
one sliding ring. The increased entropy of the segment 
on which the ring is sliding biases the amount of time 
monomers spend on that side. The solid line in Fig. 0] 
corresponds to w 10 (£,L) = Atl{l x ™{L - £) 176 ). Note 
that this functional form is still singular at I — > 0, al- 
though the singularity is weakened compared to the case 
without the sliding ring. This singularity is not visible in 
Fig. [3 presumably due to the short size of the simulated 
chains. 

The entropic effects can be increased even further by 
placing more sliding rings on the loops of Fig. [21 The 
presence of n sliding-rings on a loop enhances the number 
of configurations to 
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FIG. 4: Probability distribution of the size £ of a loop with one 
sliding ring, in competition with a simple loop. The dots are 
the result of a simulation with over 10 9 MC steps for L = 50. 



w(t,L) 

11 




20 40 60 80 100 

FIG. 5: Probability distribution of length for two loops with 
one ring on each side (L = 100). 



The side with the larger number of sliding rings is ex- 
pected to dominate in a competition. FigureElshows the 
probability distribution when we have one ring on each 
side. The solid line in Fig. [5] corresponds to the formula, 



ll'n 



t\{L-l)\ 



{I - n)\{L - I - n)!^- 76 (i - 1) 



1.76 



(5) 
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with n — \. The dots on the figures are the result of a 
simulation of 10 9 MC steps. Similar results have been 
observed for chains with sizes L — 50, 150, and 200. 
The corresponding results with two rings on each side 
are shown in Fig. EI the solid line represents Eq. JSJ with 
n = 2. Note the dramatic difference between the shapes 
of these two figures: with one ring on each side the distri- 
bution is peaked (in fact singular) at the two extremes, 
while for n = 2, the maximum has moved to the cen- 
ter. This trend becomes even more pronounced in Fig. [3 
which illustrates the case with n — 4. Clearly, increasing 
the number of rings results in a probability distribution 
which peaks more sharply in the middle. The solid line 
in Fig. again represents Eq. (JSJ) with n = 4. The good 
matches between the probability distributions obtained 
in simulations, and Eq. (|SJ) for different numbers of rings 
n lend further credence to the validity of our method. 



FIG. 6: Probability distribution for two loops with two rings 
on each side (L — 100). 
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FIG. 7: Probability distribution for two loops with four rings 
on each side (L = 100). 
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IV. KNOTS 
A. Background 

We are now in position to apply the entropic competi- 
tion method to the more complicated problem of knots. 
Knots frequently appear in closed polymers and play a 
major role in numerous biological systems. For example, 
during transcription of DNA, a variety of "de-knotting" 
enzymes, called topoisomerases, remove knots and en- 
tanglements to allow this process to go forward [r?l Hsj . 
Understanding of the action of these enzymes has been 
improved with the help of knot theory. Knotted config- 
urations have also been observed in some proteins and 
interfere with their folding into the proper shapes and 
thus lead to diseases, such as Alzheimer's and primary 
amyloidosis 19|. While knotted configurations hamper 
the proper functioning of a number of biological and 
chemical systems, they have crucial and positive roles 
in many others. The tertiary structure of RNA is an ex- 
ample in which the topological entanglements may have 
positive influences. There has been extensive research on 
predicting and understandin g th e relation between struc- 
ture and function in RNA [2(|. The folded structure 
of some types of RNA molecules (Ribozymes) determine 
their catalytic activities which are crucial to the func- 
tioning of the cell^3|- Pseudo-knots, which are formed 
by base pairing between a loop and some region out- 
side the loop, have been found in various kinds of RNAs 
and are recognized as simple RNA folding motifs j^]. A 
more physical example is provided by permanent entan- 
glements in rubber, which restrict the number of allowed 
configurations of each segment, and thus influence the 
elasticity of rubber |2^|. 

On the experimental front, artificial knots have been 
tied on both DNA and actin filaments |2^|, and the first 
Borromean DNA rings have also been assembled[24| . Sin- 
gle molecule techniques have provided a powerful tool to 
examine a wide variety of physical properties of knotted 
polymers [25|. Despite all this progress, our knowledge 
about the typical conformations and physical properties 
of knotted polymers remains rudimentary. This is due 
to the difficulty of incorporating topological constraints 
in the analytical treatments of the statistics of polymers, 
as mathematical methods of knot detection are mainly of 
"algorithmic" nature. |2(j. This complexity has encour- 
aged the use of Monte Carlo simulations, such as the one 
reported here. 



B. Competing knots 

How does the altered topology of a closed curve with a 
knot modify the number of available configurations? To 
investigate this issue, Orlandini et a/.|(| performed grand 
canonical Monte Carlo simulations of SA polygons with a 
fixed knot type K,. In these simulations the length of the 
chain is not fixed, but the set of allowed moves is such 
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FIG. 8: Probability distribution for two loops with one knot 
on each side (L = 100). 



that the topology of the chain is preserved. Orlandini 
et al. conjectured that the number of configurations of 
knotted loops takes the same general form as Eq. QJ, 
asymptotically behaving as 



(0) 



with parameters that may depend on the knot type K,. 
Assuming this form, one can find an expression for the 
mean length (n(JC)}, a quantity that can be measured in 
the grand canonical simulations. Fitting the simulation 
data for (n(/C)) _1 , and extrapolating the results to the 
limit of infinite size, Orlandini et al. can estimate the 
parameters in Eq. |JBJ. In particular they confirm that 
the effective growth constant fj,ic is independent of the 
knot type. Assuming that this is the case, they then 
conclude that the universal power-law exponents behave 
as a(/C) = a(0) + Nf, where refers to a simple loop 
(unknot), while Nf is the number of prime factors in the 
knot type /C. Such a conclusion has a simple and elegant 
interpretation: each prime factor of the knot becomes a 
tight element that incorporates an asymptotically small 
fraction of the monomers. The tight factors can occupy 
any position along the remaining large loop, each increas- 
ing the number of configurations by a factor of roughly 
the size of the loop, much like the sliding rings discussed 
in the previous section. 

If a prime knot, such as a trefoil, increases the number 
of configurations of the closed polymer by the same factor 
as a sliding ring, the two should behave similarly in the 
arena of "entropic competition." Indeed a knotted loop 
added to one side of Fig. Epulis the entire chain on its side 
as in Fig.0] However, when we place one trefoil knot on 
each of the loops, the resulting segment distribution, as 
plotted in Fig. [H] is qualitatively different from the curve 
found when a sliding ring appears on each side (Fig. EJ. 
The distribution for the competing knots is peaked at the 
center, much like cases with more than one sliding ring 
on each side. The solid line in Fig. [5] in fact represents 
Eq. JSJ with n — 4. We do not suggest that the trefoil 
is asymptotically similar to four rings, but that this is 
a good effective description of a trefoil knot with a few 
hundred monomers. To obtain more insight on finite- 
size effects, we compare the probability distributions of 
knotted loops with different sizes ranging from L = 50 to 
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FIG. 9: A logarithmic plot of probability density for two loops 
with one knot on each side for N = 50 (square), 100 (dashed 
line), 150 (dots), and 200 (solid line). The graphs were verti- 
cally displaced to coincide at £/L = 0.5. As L increases, the 
distributions become wider. 



L = 200. 

Figure [5] shows the logarithmic plot of rescaled prob- 
ability densities for N = 50, 100, 150, and 200, as a 
function of £/L. The y coordinates are shifted so that 
the maxima of distributions for all sizes coincide. We 
note that the maximum becomes flatter upon increas- 
ing L. The lack of data collapse is a clear indication of 
finite-size effects. To check this, we performed similar 
simulations with four rings on each side with chains of 
lengths L = 50 through L = 200. In contrast to knots 
(Fig.|5J), we observed that all rescaled curves collapsed in 
this case. 



V. KNOTS VERSUS RINGS 

In order to compare the relative "strengths" of knots 
and rings in increasing the number of configurations, we 
performed several simulations in which we pitted a tre- 
foil knot against different numbers of rings, as depicted 
in Fig.HHI While the hypothesis of a tight knot |1 111 li^ 
suggests that the trefoil should be well matched to a sin- 
gle ring, in our simulations with chains of 100 monomers, 
we find that around 6 rings are necessary to prevent the 
chain from being pulled completely to the side with the 
trefoil. Part of this effect is no doubt due to the size of 
the knot: even in its most compact form the trefoil knot 
involves around 14 monomers, while the sliding ring is 
assumed to occupy only one. In actuality, even a 'tight 
knot' will most likely be considerably bigger than the 
minimal size J _with typical sizes that grow with the length 
of the chain [23. To test for these effects we performed 
some further studies as reported in this section. 

To take account of the minimal size of the trefoil knot, 
we also added a similar constraint to the side with sliding 
rings. In these simulations the monomers were prevented 
to move from the ring side to the knot side with a prob- 
ability of l/(£ — 13). This means that when £ = 14, 
the monomer is strictly forbidden to pass through the 
hole from the ring side to the knot side; thus, at least 14 
monomers remain on each side throughout the simula- 
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FIG. 10: Schematic "entropic competition" between a trefoil 
knot and 6 sliding rings. 




FIG. 11: Two loops are separated by the thick oval shaped 
hole. The dotted ovals are slip-links. 



tion. Even in this case, we observe that the trefoil "wins" 
the competition against a single sliding ring. Once more, 
we found that at least four rings are need for the two 
sides of the chain to exert equal amount of force on each 
other. Thus the minimal size of the knot is not a crucial 
issue. 

We certainly do not expect the trefoil to be maximally 
tight. The number of monomers participating in the knot 
region is itself a not-so- well-defined and fluctuating quan- 
tity. A better model than the sliding ring, which allows 
for this possibility, is the slip-linked para-knot depicted 
in Fig. 1111 The thick oval ring in this figure represents 
the 'worm-hole' separating the two segments, while the 
dotted ovals are slip-links which keep two points on each 
loop close to each other (creating a figure-8 structure on 
each segment). The effect of each slip- link is similar to 
the hole in that they create two loops; however, there 
are self-avoiding interactions among all the monomers of 
a loop made into a figure-8 by the slip-link. The number 
of configurations of a figure-8 structure with self-avoiding 
constraints for 1 -C £ -C L scales as £~ C (L — £)~ c with 
c = 2.88 in d = 2 and c = 2.24 in d = 3 0. Al- 
though this formula is strictly valid only for large £ and 
L, a recent experiment on a figure-8 chain on a vibrat- 
ing plate in 2-dimensions is in good agreement with this 
formulad^lp- 

Extending the asymptotic formula for the figure 8 to 
all separations, we can obtain an analytic form for the 
number of configurations of two figure 8's in competi- 
tion. Assuming a minimal size s m i n for each segment, 
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FIG. 12: Probability distribution for two loops with one slip- 
link on each side (figure-8), from Eq. J7J with L = 100 and 
Smin = 14 (a) and s min = 30 (b). 

the number of configurations is 

Each integrand represents the probability density of a 
loop with its associated para-knot. Since there exists self- 
avoiding interaction between each loop and its para-knot, 
we set c = 2.88 and c = 2.24 for two and three dimension 
respectively. FigurelT^k is a graph of Eq. Q) with s mi „ = 
14 and c = 2.24. The distribution shown in Fig. 112b is 
qualitatively similar to the one in Fig. [SJ for two loops 
with one sliding-ring on each side, in that for both cases, 
the minimum of the distributions is in the middle and 
their maxima are on the sides. If we increase the minimal 
size from s m i n = 14 to s m i n = 30, the maximum moves 
to the middle as depicted in Fig. 112b . 

VI. CONCLUSIONS 

An important quantity in polymer physics is the num- 
ber of configurations that a chain of length I can take. As 
noted in the introduction, this quantity has the asymp- 
totic form of [i e £~ c where /i depends on the microscopic 
features of a chain while the exponent c is universal. In 
this paper we focus on obtaining the universal exponent, 
while bypassing the parameter fx. To this end, we employ 
a hole (or buckle) to divide a polymer into two segments 
(see Figs, n El i and then allow the two parts to exchange 
monomers and to "entropically" compete. The resulting 
histograms allow us to calculate the entropic exponents 
with a variety of boundary conditions and/or topological 
constraints. 



The scaling form in Eq. [5] is only valid asymptoti- 
cally, and there are in general corrections to scaling which 
complicate the extraction of entropic exponents for short 
polymers. In the application of entropic competition to 
knots, we find that finite size effects are so pervasive that 
we cannot extract reliable universal information pertain- 
ing to knots for (L < 300) . (We find an effective exponent 
of around 4 for a trefoil knot over the simulated lengths.) 
Indeed, simulations performed by Farago et a/.|28( in- 
dicate that while the ratio of the knot size N/c to the 
total length of the chain decreases, the effective knot size 
grows with the size of the polymer, suggesting that simu- 
lations on very large chains are needed to avoid the finite 
size effect. An important question which remains is at 
what size of a molecule the universal predictions of the 
asymptotic theories are expected to hold. 

The primary goal of this paper is to present the method 
of "entropic competition" for calculation of entropic ex- 
ponents of polymers with different boundary conditions 
and under different topological constraints. Our pro- 
cedure for extracting exponents is in principle simple, 
yet capable of producing accurate results. An intriguing 
application of the method to topologically constrained 
chains, such as knotted polymers, is also attempted in 
this paper. This example also illustrates the limitations 
of the technique, in that we find that finite-size effects 
are quite important to our simulated knotted chains with 
L < 300 monomers. 

It's noteworthy to mention that each simulation re- 
ported in this paper was obtained in a maximum of two 
week's CPU time on a desk-top computer. One can ap- 
ply "entropic competition" method to longer chains and 
calculate "entropic exponents" with higher accuracy on 
a large cluster of computers with longer amounts of CPU 
time. 
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